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We study a model of a 2D ultracold atomic gas subject to an “optical flux lattice”: a laser 
configuration where Raman-dressed atoms experience a strong artificial magnetic field. This leads 
to a bandstrncture of narrow energy bands with non-zero Chern nnmbers. We consider the case 
of two-level (spin-1/2) fermionic atoms in this lattice, interacting via a repulsive s-wave contact 
interaction. Atoms restricted to the lowest band are described by an effective model of spinless 
fermions with interactions that couple states in a momentum-dependent manner across the Brillouin 
zone; a consequence of the Raman dressing of the two spin states. We present the results of detailed 
exact diagonalization studies of the many-body states for a range of filling factors, v. First, we 
present evidence for the existence of a phase with coupled ferromagnetic-nematic ordering, which 
was previously suggested by a mean-field analysis. Second, we present evidence indicating the 
presence of a Laughlin-like fractional quantum Hall state occurring at filling factor v = 1/3. Finally, 
we observe a charge density wave state at iz = 1/2, which we are able to cleanly distinguish from the 
Laughlin-like state by its translational symmetry breaking and relatively small participation ratio. 

Subject Areas: Condensed Matter Physics — Strongly Correlated Materials; Atomic and Molecular 
Physics. 


There is significant interest currently in experimen¬ 
tal realizations of both topological phases of matter and 
itinerant ferromagnetism in ultracold atomic gases; 
notwithstanding that, to date, topological phases have 
not yet been observed,while ferromagnetic behaviour 
has proved challenging to detect.In this paper we 
present evidence demonstrating the promise of “optical 
flux lattices” (OFLs) — optical lattices that generate ar- 
tihcial magnetic helds via coherent Raman coupling of 
the internal states — to advance both aims.^^’^^ We do 
so by presenting the first results from exact diagonal¬ 
ization calculations of strongly interacting fermions in 
optical flux lattices. (Previous studies have considered 
bosonic atoms. 

The key feature that leads to the usefulness of OFLs in 
stabilizing both ferromagnetic phases and topologically- 
ordered phases is that, throughout a large region of their 
parameter space, they exhibit narrow, almost flat bands 
with non-trivial band topology (non-zero Chern number), 
due to the presence of a net non-zero artihcial magnetic 
field. 

From the perspective of ferromagnetic ordering, the 
band flattening leads to an enhanced density of states 
without exponentially localized Wannier orbitals,al¬ 
lowing the Stoner instability to a ferromagnetic phase to 
occur at a value of interaction strength that is much re¬ 
duced from the value required in free space. This has 
been demonstrated within a mean-held treatment of the 
interactions,^® which also shows that, owing to the Ra¬ 
man coupling of the spin states, the ordered phase has 
combined ferromagnetic-nematic order, with broken ro¬ 
tational symmetry of both the spin and orbital motion.^® 
Given that this phase appears in the strongly interacting 
regime, it is important to question the extent to which 
mean-held theory accurately describes the system. Our 


exact diagonalization studies provide clear evidence for 
the presence of ferromagnetism and the associated ro¬ 
tational symmetry breaking. Indeed, we hnd that the 
ferromagnetic-nematic phase appears at an even smaller 
interaction strength than in mean-held theory. The in¬ 
teraction strength is critical to practical implementations 
because it greatly inhuences the rate of three-body loss 
processes, which amongst other challenges, impede ex¬ 
perimental progress.That the interaction strength 
is much reduced from its free-space value by the OFL 
can very signihcantly reduce these loss processes. 

From the perspective of topological phases, the nearly 
hat and topological bands of the OFL make them ideal 
candidates in the search for topological fermionic phases 
analogous to the fractional quantum Hall (FQH) effect 
of electrons. Here we shall provide compelling evidence 
demonstrating that indeed there is a strong FQH-like 
phase analogous to a Laughlin state at hlling factor v = 
1/3 for interacting fermions in the lowest band. There 
are a number of existing theoretical proposals to engi¬ 
neer fermionic FQH-like phases in band structures with 
non-trivial topology, for instance using a combination 
of precisely tuned nearest-neighbour and next-nearest- 
neighbour tunnelling and either on-site interactions, 
or dipolar interactions.^® Such phases are often termed 
fractional Chern insulators (FCI). The presence of a 
Laughlin-like state in an OFL is intriguing in this con¬ 
text because the interactions have a rather different form 
from what has been considered previously. Although 
they arise from a 2-body s-wave contact interactions, 
once projected to the lowest band they acquire a non¬ 
local character with an unusual k-dependent structure 
(a feature that originates from the optical dressing as we 
shall describe in more detail momentarily). Nevertheless, 
we find that the OFL’s Laughlin-like state is surprisingly 
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robust. 

In addition to the aforementioned ferromagnetic- 
nematic phase and the FQH-like phase, we observe a 
highly incompressible state at v = 1/2. Unlike at i' = 
1/3, we show that this state breaks both translational 
and rotational symmetry, which allows us to clearly dis¬ 
tinguish it as a charge density wave (CDW). 

Model 

An optical lattice is formed by the interference pattern 
of linearly polarized in-plane laser beams with common 
wave vector, n. We study a hexagonal lattice formed 
by three such beams at orientations differing by 27r/3, 
which are described by wave vectors Ki = —k/2{^/3, 1), 
K 2 = k/ 2(-\/3, —1) and K 3 = k( 0, 1). In order to engi¬ 
neer a spin-like degree of freedom one exploits the fact 
that atoms possess hyperfine atomic levels (which can be 
split by the Zeeman effect). For simplicity we focus on a 
2-level system with angular momentum F = 1/2. Prac¬ 
tical implementations might have higher values of F, but 
the key features at F = 1/2 should remain.^® The 2-level 
system is coherently Raman coupled, or dressed, by a 
fourth, circularly polarized light beam oriented perpen¬ 
dicular to the 2D plane of the lattice. This configuration 
leads to a 2-component Hamiltonian (in the notation of 
Ref. 18) 

^0 = + Ko(r)i + d- • B(r). (1) 

The optical lattice is described by a scalar potential 

Kc = Vo [3cos^(0) — 1 ] '^2 cos(k'- • r), (2) 
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and the Raman coupling is described by a scalar product 
of a vector of Pauli matrices a- with a local vector field 
B, whose components are 

Bz = VSVo sin^(d)^ sin(K' • r), (3) 

3 

Boo + iBy = eVo cos( 6 ») ^ (4) 

3 

where /«) = H 2 — >^ 3 , K 2 = K 3 — Ki and = Ki — 

The tunable model parameters are the lattice depth, Vq, 
the linear polarization angle 9 of the beams forming the 
optical lattice, and the amplitude ratio e between the 
in-pane beams and the fourth, perpendicular beam. To 
reduce the number of free parameters, and for direct com¬ 
parison with the existing mean-field theory results, we fix 
e = 0.4 and 9 = 0.3.^® We shall, however, consider the 
dependence of the properties on lattice depth Vq, as com¬ 
pared to the recoil energy, En = /2M. 

The spectrum of Eq. 1 exhibits a set of energy bands 
with topological character. In particular, for any non¬ 
zero Vo! the lowest energy band has unit Chern num¬ 
ber. It is analogous to the lowest Landau level in a sys¬ 
tem in which there is one flux quantum per real space 


unit cell.^"^ Thus, the “quantum Hall” filling factor, v — 
defined as the ratio of 2D particle density to magnetic 
flux density — is equivalent to the band filling factor, 
i.e. the number of atoms per unit cell. The dispersion of 
this lowest energy band e(k) varies with Vo/En, having 
a minimum at around Vo/En = 2 .^^ 

For our numerical studies, we construct the energy 
eigenstates of Eq. 1 by proposing a set of Bloch wave 
functions on a Li x L 2 k-space grid {K} = k[Ii/Li — 
K 2 I 2 /L 2 for Ij = 0 ,... Lj — 1 with j = 1 , 2 , and a set of 
reciprocal lattice vectors {G}: 

Y. (5) 

g6{G} 

where n is a band index, ct is a spin index and k G {K}. 
The Raman coupling induces a k-dependent mixture of 
spin-up and spin-down character, € [— 1 / 2 , 1 / 2 ], in 
the single-particle states as a function of the amplitudes 
^k CT- the hexagonal Brillouin zone of the lowest band, 
single-particle states close to the K/K' points obtain a 
largely spin-up/spin-down character. The interaction- 
driven ferromagnetic-nematic phase is, therefore, asso¬ 
ciated with a breaking of the original 6 -fold symmetry 
of the Brillouin zone down to a 3-fold symmetry where 
states close to K jK' are unevenly populated. 

The simplest interactions that can be introduced 
into this model are s-wave contact interactions be¬ 
tween spin-up and spin-down. We represent these by 
a two-dimensional delta-function interaction potential, 
of strength /M, setting the dimensionless interaction 
strength g. We consider repulsive interactions (positive 
g). In practical implementations, the magnitude and 
sign of these interactions can be tuned using a Eeshbach 
resonance.^® The projection of that contact interaction to 
the lowest band produces an effective interaction between 
spinless fermions (again, owing to the Raman coupling). 
The effective Hamiltonian describing the lowest band is 

Bn^ = ^ T 2 ^ ^ hk^k2k3k4C^ki C^k2^k3Ck4 , 

k kik2k3k4 

( 6 ) 

with 

Vkik2k3k4 = ^ [ c^^r5I'^ki^(i')</'k2a(r)</k3a(i')</k4<T(r), 

( 7 ) 

where cr denotes a spin flip, and we have suppressed the 
band index n since we have projected to the lowest band. 

We shall present detailed analyses of the ground states 
and low-energy spectra of iLib for a range of filling fac¬ 
tors, lattice depths and interaction strengths. We study 
N particles on the simulation grid {K}. One is free to 
change the offset Ai or A 2 between the grid {K} and the 
underlying Brillouin zone along directions and —K 2 
respectively. (This freedom is equivalent to flux inser¬ 
tion in the two cycles of the torus formed by the periodic 
boundary conditions.) By varying these offsets we can, 
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for example, check the extent to which our results are 
influenced by hnite-size effects. 


Ferromagnetic—nematic order 

First, we shall identify the presence of a ferromagnetic- 
nematic phase in this model at low Filing factors. That 
two forms of ordering appear in this phase is a conse¬ 
quence of the spin-orbit coupling inherent in the Raman 
coupling. 

To quantify the appearance of ferromagnetic-nematic 
order, we measure the normalized density-density corre¬ 
lation function weighted by the spin character Sk at each 
point in fc-space, given by 

X! (‘S'k'c'''k'Ck'S'kc'fkCk) • ( 8 ) 

k,k' 

For complete ferromagnetic order, with ^k = -1-1/2 or 
S'k = —1/2, one would have S'kO^kCk = N/2 and ps = 

l. Note, however, that the spin S'k has some dependence 
on the lattice depth parameter, Vq, and |Sk| < 1/2. 

In order to diagnose the presence of (or lack of) the 
associated nematic order, we study rotational symmetry 
breaking in fc-space. We calculate a rotational density- 
density function that measures commensurablility under 
changes in angular momentum m: 

^ (^^k'^AokCkCk) , (9) 

k,Q! 

where Ra are operations to rotate the k-states by angle 
a. In order to elicit signatures indicating breaking of 
the original 6 -fold symmetry we select a = titt/S and 

m, n G {0-5}. Pr(to) is only defined when the boundary 
conditions of the simulation grid themselves obey 6 -fold 
rotational symmetry. The signature of fc-space rotational 
symmetry breaking is that fluctuations 

^ PrM-(PrM) 

dp^{m) = - -——y - (10 

in /Or(to) satisfy |( 5 / 9 R(m)| ~ 1 for some m. Note that we 
choose a convention where the mean value (pR(m)) does 
not include Pr( 0) but does include Pr(1) to Pr(5). 

In Fig. 1 we plot a map of the spin-weighted density- 
density function (Fig. la) and a map of the fluctu¬ 
ation in the rotational correlation function 5p^{m = 3) 
(Fig. lb) for the ground state at Ailing factor v = 1/A for 
a range of interaction strength g and lattice depth Vq. 
The plots are for particular values of the offsets Ai and 
A 2 that allow {K} to have 6 -fold rotational symmetry. 
We have also checked that the values of ps are largely 
independent of the boundary conditions. Non-zero ps in¬ 
dicates a magnetic ground state and 6p^{3) ^ 1 indicates 
that the 6 -fold rotational symmetry is broken to 3-fold. 



FIG. 1. (color online). Map of the ground state character 
of Rib (Eq. 6) as a function of lattice depth Vo and interac¬ 
tion strength gatj/ = l/4[N = 9ona6x6 simulation 
grid]: (a) Spin-weighted density-density function ps [defined 
in the text], indicating a (partially) spin-polarized ground 
state where 1 > ps > 0. (b) Fluctuation (5pr(3) of Pr(3) 
from its mean value, as defined in Eq. 10, indicating breaking 
of 6-fold to 3-fold rotational symmetry where (5 pr(3) ~ 1. 


From these results we conclude that evidence from 
these small systems supports the mean-field theory pre¬ 
diction of interaction-induced ferromagnetic-nematic or¬ 
der. Our results differ from mean-field theory in that 
they typically predict a partially polarized, rather than 
fully polarized ferromagnet. Also, the value of g nec¬ 
essary to see some degree of ferromagnetism is signifi¬ 
cantly smaller than the mean-field theory prediction (see 
Ref. 18). Indeed, from Fig. 1 it would appear that rela¬ 
tively weak interactions and small Vq s-re sufficient for fer¬ 
romagnetism to occur. In particular, for non-zero VojEji, 
the ferromagnetic-nematic phase appears in regimes of 
interaction g that are as much as ten times smaller than 
the value p = 27r of the Stoner instability for 2D itinerant 
fermions in free space. This observation is very encour¬ 
aging from an experimental perspective, as it indicates 
that ferromagnetic ordering can be found in regimes of 
interactions far from the Feshbach resonance where three- 
body losses are much suppressed. 


Strongly correlated phases 

Next we turn our attention to the identification of 
strongly correlated phases in the OFF. In Fig. 2 we plot 
the inverse compressibility /3~^{h') for the ground state of 
the lowest OFF band as a function of Ailing factor ^ when 
the model is in its narrow band regime. This quantity is 
given by estimating the second derivative with respect to 
particle number N of the ground state energy 

= N/[L,L2]) = ^ [E^+^ + Eg-^ - 2Eg] . 

( 11 ) 
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FIG. 2. (color online). Inverse compressibility P~^[u) [de¬ 
fined in the text] for the ground state of Hi^ (Eq. 6) as a 
function of filling factor v of the lowest OFL band, for dif¬ 
ferent values of interaction strength g in the narrow band 
regime. Calculation for a 5 x 6 grid with VoIEr = 2. Highly 
incompressible states occur at :/ = 1/3 and = 1/2 once 
interactions are present. 


The most prominent features in Fig. 2 are peaks in the 
inverse compressibility at filling factors i/ = 1 /2 and v = 
1/3 that occur only once interactions are introduced {g > 
0 ). 

There are a number of possibilities that would be con¬ 
sistent with these inverse compressibility signatures: For 
instance a FQH-like state, a CDW, or a striped phase 
of the so-called “thin-torus pattern” (that is, stripes of 
a FQH-like ground state that would occur in the long, 
thin torus limit of the model). To assist in distinguish¬ 
ing between these different candidates we calculate the 
“fractional participation ratio”, 




( 12 ) 


where the sum is over D possible Fock states i = 1,... ,D 
with amplitudes P provides an estimate of the 

fraction of Fock states that significantly contribute to a 
many-body wave function. In Fig. 3 we show the frac¬ 
tional participation ratio for the ground state at filling 
factor z/ = 1/3. For the set of model parameters used in 
this calculation, the OFL exhibits a narrow band when 
y^/ER « 1.75-3.0. We observe that P is close to 0.33 
in the narrow band regime, for sufficiently strong inter¬ 
actions. This large value indicates that the ground state 
at = 1/3 is strongly correlated and cannot be well rep¬ 
resented by a small number of Slater determinants. On 
the other hand, we find that the fractional participation 
ratio is much smaller for the ground states v = 1/2 
and V = 1/4, taking values on the order of P = 0.01 
(for a 4 X 6 system with N = 12) and P = 0.001 (for a 
6 x 6 system with N = 9), respectively. Therefore neither 
are strongly correlated. The result that ly = 1/4 is not 
a strongly correlated phase is consistent with our con¬ 
clusion (above) that the ferromagnetic-nematic ground 
state that we have identified here using exact methods is 
close to its mean-field description. 



FIG. 3. (color online). Fractional participation ratio P 
(defined in the text), estimating the fraction of Fock states 
contributing to the ground state of Hn, (Eq. 6), for a range of 
lattice depth and interaction strength for = 10 on a 5 x 6 
grid (z/ = 1/3). 


FQH-like state at u = 1/3 


We shall now characterize the z/ = 1/3 state in some 
more detail. First we shall investigate the effect of inser¬ 
tion of magnetic flux quanta through one cycle of the 
torus. The addition of a single flux quantum, = 1, 
is simulated by changing the boundary condition A 2 by 
precisely one lattice spacing, k^IEi (or similarly for the 
/«/ direction). This procedure also shifts the spectrum 
with the same periodicity due to the underlying OFL 
band structure. In Fig. 4 we plot the energy spectrum of 
the 5x6 system as a function of N^. Owing to the pe¬ 
riodicity of the OFL in real space, the eigenstates of Pl^^ 
can be subdivided into different linear momentum sec¬ 
tors ] along the two grid directions in fc-space. 

We observe a 3-fold degenerate ground state manifold in 
the sectors expected for the Laughlin state. Once the 
underlying variation is removed—by subtracting off the 
mean energy of the ground state manifold Eo{N ^)—we 
clearly observe the characteristic spectral flow pattern of 
the Laughlin state, whereby the ground states only return 
to their original ordering after the insertion of precisely 
3 flux quanta (Fig. 4 inset). 

In Fig. 5a we plot the finite-size scaling of the many- 
body energy gap Amb for the z^ = 1/3 ground state. In 
Fig. 5b we plot the finite-size scaling of the ratio of the 
many-body energy gap Amb to the width d of the ground 
state manifold. These values are estimated by the min¬ 
imum possible Amb and the maximum possible S for a 
range of offsets Ai and A 2 . With the evidence avail¬ 
able we conclude that the energy gap remains non-zero 
in the thermodynamic limit. For the Ly = Q data set we 
observe that Amb/^ becomes negligible for larger N. In 
the “thin torus” case {Ly = 3), however, we observe that 
the ratio Amb/iJ tends to a non-negligible value. This 
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FIG. 4. (color online). Energy spectrum of flib (Eq. 6) for 
filling factor v — 1/3 of the lowest OFL band [A^ = 10 on a 
5x6 grid] with VolEa = 2.0 and g 1271 = 1 as a function the 
number, N^, of flux quanta inserted along the K 2 direction 
= 1 corresponds to A 2 = K 2 /^ 2 ). Inset: Magnification 
of the ground state manifold over the same region with the 
mean ground state energy, Eo{N^), subtracted off. The 3 
ground states exhibit spectral flow whereby they return to 
their original ordering only after 3 flux quanta are inserted. 


A„b /Efi (a) 


5/A„b (b) 



FIG. 6. (color online). Energy spectrum of Hib (Eq. 6) 
for A^ = 9ona5x6 grid with VoIEr = 2.0 and g/27v = 1 
as a function of linear momentum sector The 

counting of quasi-hole states (energy levels below the gap) 
per sector precisely matches the analytic result given e.g. in 
Ref. 22. 
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FIG. 7. (color online). Ground state character of Hu, 
(Eq. 6) for filling factor v = 1/3 of the lowest OFL band 
[A^ = 12 on a 6 X 6 grid]: (a) Representative fluctuations 
Spt{ci) in the static structure factor pT(q) [defined in the 
text]. ]5pT(q)| ~ 1/A^ here indicates translational invariance; 
(b) Representative fluctuations SpR(m) in the rotational cor¬ 
relation function puirn) [defined in the text]. [5pR(m)[ ~ 1/N 
here indicates 6-fold rotational invariance. 


FIG. 5. (color online). Thermodynamic extrapolation of 
the many-body energy gap Amb (defined in the text) of 7/ib 
(Eq. 6) at filling factor = 1/3 for different aspect ratios 
of the simulation grid and for g/27i = 1. Dashed lines show 
linear extrapolations, where helpful: (a) Many-body energy 
gap, Amb. (b) Ratio of the energy spread, 5 (defined in the 
text) to the energy gap Amb . 


behaviour is similar to that reported for FCI models. 

The spectrum of quasi-hole states can be obtained by 
removing particles from the 1/3 state (for instance, N = 
9 in a 5 X 6 system, see Fig. 6) and we find that the 
counting of quasi-hole states per linear momentum sector 
precisely matches the analytic formula given in Ref. 22. 

To rule out the possibility that the v = Ij?) state could 
be a CDW or striped phase, we look for the presence or 
absence of translational symmetry breaking by comput¬ 


ing the static structure factor 

PT(q) = ^ 'y ^ (c^k'-qCk'C^k-hqCk) ) (13) 

^ k,k' 

where q G {K}. As with the rotational correlation func¬ 
tion, we can identify translational symmetry breaking by 
analysing the order of magnitude of fluctuations from the 
mean value 

IT PT(q) - (pT(q)) ,,,, 

= (pM) ' 

where again we choose the convention to disregard Pt (0) 
from the calculation of the mean. CDW or striped states 
break translational symmetry and will satisfy |<5pT(fl)| 

1 for some q, but FQH-like states do not. FQH-like states 
should also be rotationally invariant, when possible on a 
lattice, owing to their uniform density. 

In Fig. 7 we plot the fluctuations in the static structure 
factor (Spt) and also the rotational correlation function 























































6 


8.0 




6.0 


g 

tq 

I 


4.0 


2.0 


------ **±SX*********3 

.***♦**•**= 



****♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦-♦*♦** 


0.1 0.2 0.3 0.4 


0.5 


0.6 0.7 0.8 0.9 1.0 


FIG. 8. (color online). Energy spectrum of Hu, (Eq. 6) for 
filling factor v = 1/2 oi the lowest OFL band {N = 12 on a 
4x6 grid) with Vo/Ea = 2.0 and g/2n = 1, as a function the 
number, N^, of flux quanta inserted along the direction 
= 1 corresponds to A 2 = n^/L-^). For clarity, for each 
the minimum energy, Eo(A^,^), has been subtracted off. 
No simple spectral flow is observed. 


the torus. No simple spectral flow is observed and in¬ 
stead we see dispersion and mixing of a large number of 
low-lying energy levels. In Fig. 9a we plot representa¬ 
tive fluctuations in the static structure factor {5p^) as¬ 
sociated with these low lying states. Clear peaks are 
seen in Sp^ corresponding to momentum sectors satisfy¬ 
ing mod(gj, Lj/2) = 0 for j = 1, 2. This character is also 
observed in the ground state and low lying states of the 
V = 1/2 state in 4 x 4, 4 x 8 and 5x6 systems. It indi¬ 
cates the formation of a CDW ground state with a unit 
cell that is doubled in size along both axes. In Fig. 9b 
we plot representative fluctuations in the rotational cor¬ 
relation function ((5 pr) associated with these low lying 
states. We observe that they retain the 6-fold rotational 
symmetry in fc-space. Our conclusion here is consistent 
with the above indication from the fractional participa¬ 
tion ratio that ;/ = 1/2 is not strongly correlated. 


Conclusion 
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FIG. 9. (color online). Gharacter of the low lying states 
of Hu (Eq. 6) for filling factor v = 1/2 of the lowest OFL 
band with Vo/Er = 2.0 and g/2'iT = 1: (a) Representative 
fluctuations 5pT(q) in the static structure factor pT(q) [de¬ 
fined in the text] for the lowest lying states {N = 12 on a 
4x6 grid). Peaks with |(lpT(q)| ~ 1 here indicate trans¬ 
lational symmetry breaking; (b) Representative fluctuations 
Spuijn) in rotational correlation function puijn) [defined in 
the text] for the lowest lying states {N = 8 on a 4 x 4 grid). 
]5pR(m)] ~ l/N here indicates 6-fold rotational invariance. 


(i5pr) for the ground state at = 1/3. We clearly observe 
that it breaks neither translational nor (6-fold) rotational 
symmetry in A:-space. The structure of the fluctuations 
is associated with the properties of the quantum Hall 
liquid (see e.g. Ref. 31 and 32 for a discussion of the sub¬ 
structure of the static structure factor for a FQH state). 
Our study is restricted to very small systems and, as 
such, we cannot resolve that sub-structure in detail. 


CDW state at v = 1/2 

Finally, we focus on the incompressible phase that ap¬ 
pears at filling factor v = 1/2. In Fig. 8 we plot the 
energy spectrum of the 4x6 system as a function of 
the number of flux, inserted through the K 2 cycle of 


We have presented a significant body of evidence 
that together supports the conclusions that interacting 
fermions in an OFL can exhibit: 1). A (partially polar¬ 
ized) ferromagnetic-nematic phase whose onset occurs at 
a much reduced interaction strength compared to the free 
space Stoner criterion; 2). A Laughlin-like state at 1/3 
filling; 3). A CDW state at z/ = 1/2 in the narrow band 
regime. We have employed a number of simple corre¬ 
lation functions with which to measure the symmetries 
of the model, and in particular to identify when these 
symmetries are broken, and how strongly. These cor¬ 
relation functions allow us to distinguish a very strong 
FQH-like state from CDW and striped states. Further 
we have demonstrated that the many-body participation 
ratio measure can also be used to support such distinc¬ 
tions. The specific evidence we have gathered is for a 
particular set of OFL parameters (fixed e and 0), how¬ 
ever there is nothing special about our choice of these 
parameters and the narrow band regime pervades a large 
region of the OFL parameter space. Our results demon¬ 
strate that an OFL would be a fruitful approach with 
which to enhance ongoing efforts towards observing both 
ferromagnetic behaviour and more exotic phases of mat¬ 
ter in cold atomic gases. 
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